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ABSTRACT 


Interfacial  tensions  between  immiscible  homopolymers  are  measured 
using  an  automated  pendant  drop  apparatus,  which  utilizes  Video  Digital 
Image  Processing  techniques.  A  recently  developed  robust  shape  analysis 
algorithm  is  used  to  analyze  the  experimental  drop  profiles.  The  data  show 
the  effect  of  temperature  and  number  average  molecular  weight  (Mn)  on  the 
interfacial  tension  for  the  immiscible  blends  polystyrene-poly(methyl 
methacrylate),  polybutadiene-poly(dimethyl  siloxane)  and  polystyrene-hydro¬ 
genated  1,2  polybutadiene.  Interfacial  tension  decreases  linearly  with 
temperature  and  increases  with  molecular  weight.  The  data  for  all  three 
systems  can  be  approximated  by  an  empirical^M^^V^lationshlp. 

The  interfacial  tension  data  for  the  latter  blend  system  are  compared 

with  thermodynamic  theories  of  polymeric  interfaces.  A  square  gradient 

theory  approach,  in  conjunction  with  the  Flory-Huggins  expression  for  the 

free  energy  of  mixing,  predicts  a  magnitude  and  temperature  dependence  of 

interfacial  tension  which  are  In  reasonable  agreement  with  experimental 

data,  and  does  predict  a  molecular  weight  dependence  which  roughly  follows 
-2/3 

the  empirical  M  '  relationship.  " 


* 

I 

[ 


INTRODUCTION 


The  structure  and  thermodynamic  state  of  polymeric  interfaces  are 
important  features  in  many  polymeric  materials  of  current  technological 
interest.  This  is  especially  true  for  multiconstituent  systems  such  as 
immiscible  polymer  blends  or  microphase  separated  copolymers,  where  inter¬ 
face  structure  can  affect  greatly  the  mechanical  properties^.  Interfacial 
tension  is  important  due  to  its  influence  on  the  morphology  of  multiphase 
polymers.  A  number  of  experimental  investigations,  for  example,  have  shown 
that  the  phase  structure  (e.g.  dispersed  particle  size)  in  incompatible 
polymer  blends  is  directly  proportional  to  the  interfacial  tension.  q 

Fundamentally,  interfacial  tension  is  a  thermodynamic  property  of  the 
system  which  may  be  calculated  directly  from  statistical  thermodynamic 
theories.  Experimental  measurement  of  interfacial  tensions  is  therefore  a 
straightforward  means  for  evaluating  the  validity  of  these  theories. 

The  existent  data  pertaining  to  the  interfacial  properties  of  multi¬ 
phase  polymers  have  been  well  summarized  in  several  reviews,®’®^  and  a 

Q 

monograph  by  Wu  .  Interfacial  tension  generally  decreases  linearly  with 

_2 

temperature  with  a  temperature  coefficient  in  the  range  of  10  dyn/cm°C. 

An  increase  in  molecular  weight  leads  to  an  increase  in  the  interfacial 

tensions;  however,  there  are  only  a  few  studies  which  have  examined  this 

9  10 

dependence  in  any  detail.  Experimental  interfacial  tensions  between 

n-alkanes  and  a  per-fluoroalkane  (C^  5^27^’  P°ly(dimethyl  siloxane)  (PDMS) 

and  C12  5^27  or  ^18’  and  a^kanes  polyethylene  glycol)  all  exhibit 
-2/3 

an  apparent  Mn  dependence  on  molecular  weight  (where  Mn  is  the  number 
average  molecular  weight),  similar  to  what  has  been  found  for  homopolymer 
surface  tension®. 


-  2  - 


This  similarity  is  predicted  by  several  empirical  theories  which 

n  12 

relate  interfacial  tension  to  the  pure  component  surface  tensions.  ’  A 

1 3-20 

number  of  thermodynamic  theories  have  appeared  that  predict  the 

interfacial  tensions  of  polymer  blends.  The  applicability  of  these  theories 

5  7  8  21  22 

has  been  discussed  in  detail  in  several  review  articles.  *  *  *  '  To 

date,  these  theories  have  not  been  compared  to  data  with  systematic  changes 
in  molecular  weight  and  temperatures.  In  general,  the  theories  either 
assume  infinite  molecular  weight,  or  do  not  reproduce  the  observed  Mn~^3 

dependence.  Several  theoretical  treatments,  valid  near  the  critical 

1/2  19  22 

temperature,  lead  to  a  predicted  Mn  dependence  . 

In  this  paper,  the  effect  of  temperature  and  molecular  weight  on  the 

interfacial  tension  is  studied  for  three  binary  polymer  blend  systems: 

polystyrene  (PS)  with  poly(methyl  methacrylate)  (PMMA);  poly(dimethyl 

siloxane)  (POMS)  with  1,4  polybutadiene  (PBD);  and  polystyrene  (PS)  with 

hydrogenated  1,2  polybutadiene  (PBDH).  The  experimental  data  for  PS/PBDH 

blends  are  compared  to  the  predictions  of  Helfand  and  coworkers13’14  and  to 

7  ?  OA 

a  theory  based  upon  an  extension  of  the  square  gradient  theory  . 

EXPERIMENTAL 

Materials:  Trimethylsi loxy  terminated  PDMS  (Petrarch  Systems,  Inc.),  1,4 
polybutadiene  (Polysciences,  Inc.),  poly(methyl  methacrylate)  and  poly¬ 
styrene  (Polysciences,  Inc.  and  Polymer  Laboratories,  Ltd.)  were  used  as 
received.  Their  characteristics,  quoted  from  the  suppliers,  are  listed  in 
Table  1.  The  PMMA  was  prepared  by  group  transfer  polymerization  and  was 
provided  by  Or.  R.  Khanna  of  E.I.  DuPont  de  Nemours  4  Co.  The  1,2  poly¬ 
butadiene  polymer  was  anionically  synthesized  in  our  laboratory.  The 
reaction  was  carried  out  at  room  temperature  in  benzene  under  high  purity 
argon  using  a  glass  manifold  equipped  with  teflon  valves.  n-Butyllithium 
(n-BuLI)  was  the  initiator  and  dipiperidinoethane  was  used  as  a  modifier 
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at  a  molar  ratio  of  10:1  (modifier: initiator)  to  produce  homopolymers  with 
mainly  1,2-  butadiene  isomer.  A  homogeneous  catalytic  reaction  was  adapted 

nc 

from  a  method  described  by  Falk  to  saturate  the  butadiene  sequences.  The 

catalyst  used  was  a  complex  of  n-BuLi  and  the  cobalt  salt  of  2-ethyl  - 

hexanoic  acid  in  cyclohexane.  A  molar  ratio  of  Li/Co  of  2.5  was  found 

necessary  for  the  reaction  to  go  to  completion.  The  hydrogenation  reaction 

was  carried  out  at  50°C  and  50-60  psi  reactant  gas  pressure  with  stirring 

for  1-4  hours.  More  details  on  the  synthesis  procedure  are  given  by 
27 

Quan  .  The  microstructure  of  the  butadiene  sequence  and  the  extent  of  the 

hydrogenation  reaction  were  evaluated  by  proton  nuclear  magnetic  resonance 

spectroscopy.  The  Mn  and  the  molecular  weight  distributions  were  estimated 

27 

using  gel  permeation  chromatography  using  a  calibration  curve  by  Quan  . 
The  characteristics  of  the  PBDH  4080  are  given  in  Table  1. 

The  densities,  necessary  to  calculate  interfacial  tension,  were 
measured  for  PDMS,  PBD  and  PBDH  with  a  digital  density  meter  (Mettler 
Instrument  Corporation,  model  DMA45,  DMA46  and  external  cell  DM4512) 
capable  of  measuring  density  as  a  function  of  temperature  to  five  signi¬ 
ficant  digits.  The  densities  depend  approximately  linearly  on  temperature. 
The  results  of  linear  regression  are  expressed  in  the  following  empirical 
equations  which  are  used  for  the  purpose  of  interpolation: 

PDMS  770:  p  =  0 .94762-9.968. 10'4  t 
PDMS  1250:  p  =  0.96093-9. 158. 10"4  t 
PDMS  2000:  p  =  0. 98006-9. 384. 10*4  t 
PDMS  3780:  p  =  0.99773-1 .055. 10*3  t 
PBD  1000:  p  =  0.88787-5. 358. 10'4  t 
PBDH  4080:  p  =  0.88395-5.267. 10‘4  t 

3 

with  p  in  g/cm  and  t  in  °C. 
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The  density  data  of  Bender  and  Gaines  and  Fox  and  Flory  are  used  for 

30 

PS.  For  PMMA,  the  density  data  of  Somani  were  employed. 

Interfacial  Tension  Measurements 

Interfacial  tensions  were  measured  using  an  automated  pendant  drop 

31 

apparatus  which  utilizes  video  digital  image  processing  techniques  .  A 
fluid  pendant  drop  of  the  more  dense  constituent  is  formed  at  the  tip  of  a 
glass  capillary  tube  of  a  Drummond  positive  displacement  syringe.  The 
syringe  is  inserted  into  a  fluid  matrix  of  the  second  polymer  which  is 
contained  in  a  quartz  microcuvette.  The  microcuvette  is  placed  in  a 
Rame-Hart  C-2033  environmental  chamber  under  argon  atmosphere.  The  tem¬ 
perature  of  the  chamber  is  controlled  to  within  1°C. 

The  optical  system  consists  of  a  Questar  Ml  microscope  coupled  to  an 
NEC  TI-22A  CCD  video  camera.  The  optics  are  focussed  by  optimizing  the 
video  image  of  a  reticle  containing  a  finely  ruled  grid  that  is  placed  at 
the  drop  location.  The  grid  also  serves  the  in  determination  of  the  actual 
magnification  of  the  instrument. 

The  video  image  of  the  drop  is  digitized  by  a  Tecmar  Video  Van  Gogh 

frame  grabber  board  resident  within  an  IBM  XT  microcomputer.  The  resulting 
.  32 

image  is  analyzed  using  global  thresholding  and  near-neighbor  analysis  to 
extract  the  experimental  profile. 

The  equilibrium  profile  of  a  pendant  drop  is  described  by  the  equation 

33 

of  Bashforth  and  Adams 

_ 1_  +  sin<t> 

(R/a)  (x/a) 


2  +  B  (z/a)  = 


(1) 


The  shape  of  the  drop  is  specified  in  terms  of  the  x  and  z  coordinates, 
where  $  is  the  angle  between  a  tangent  to  the  drop  profile  and  the  hori¬ 
zontal  axis,  R  is  the  radius  of  curvature  at  a  point  (x,z),  and  a  Is  the 
radius  of  curvature  at  the  drop  apex.  The  shape  factor  B  is  defined  as 

B  =  f  2) 


where  y  is  the  interfacial  tension,  g  is  the  gravitational  constant,  and  Ap 
is  the  mass  density  difference  across  the  interface. 

The  experimental  profile  is  analyzed  with  a  robust  shape  comparison 
algorithm  which  is  resistant  to  outlying  points  that  may  result  during 
the  computerized  drop  profile  discrimination  procedure.  In  addition  these 
routines  minimize  the  number  of  parameters  that  must  be  searched  numeri¬ 
cally  in  the  subsequent  optimization  (i.e.  shape  comparison).  Fig.  1  shows 
the  quality  of  the  fit  obtained  for  a  PS  drop  in  a  matrix  of  PBDH. 

RESULTS  AND  DISCUSSION 

Interfacial  Tension  Data 

The  interfacial  tension  data  for  the  four  PBD/PDMS  pairs  as  a  function 

of  temperature  are  shown  in  Fig.  2.  Interfacial  tension  decreases  almost 

linearly  with  temperature  in  the  range  25°C  to  125°C.  The  temperature 

coefficients  (Table  2)  are  of  the  same  magnitude  as  those  reported  for 

34 

other  polymer  pairs  .  The  data  also  exhibit  an  increase  in  interfacial 
tension  as  the  molecular  weight  of  the  polymer  Increases.  Fig.  3  illus¬ 
trates  the  effect  of  PDMS  molecular  weight  on  interfacial  tension  (at  25°C) 
for  constant  molecular  weight  of  PBD  (M  =  982).  These  data  and  data  at 

other  temperatures  are  essentially  linear  when  plotted  against  Mn"2//3  as 

9  10 

suggested  by  Gaines  and  co-workers  ’  and  discussed  in  the  Introduction. 


The  points  for  the  lowest  molecular  weight  PDMS  specimens  (PDMS  770) 
do  not  fall  on  the  same  line.  This  may  result  from  the  ideal  phase  separa¬ 
tion  assumption  employed  to  estimate  the  density  difference  required  for 
the  interfacial  tension  calculation.  For  this  very  low  molecular  welqht, 
partial  miscibility  occurs,  even  at  the  temperature  of  25°C.  The  ideal 
density  difference  used  in  the  calculation  is  thus  too  large,  and  leads  to 
an  overestimation  of  the  interfacial  tension. 

The  interfacial  tension  data  for  the  PS/PBDH  system  appear  in  Fig.  4. 

A  linear  decrease  with  temperature  and  an  increase  with  Mp  is  also  ob¬ 
served.  The  least  squares  linear  fits  of  the  data  are  given  in  Table  3. 

The  interfacial  tension  data,  for  various  PS  molecular  weight  at  171°C  and 

a  fixed  molecular  weight  of  PBDH  also  exhibit  a  linear  relationship  when 

-2/3 

plotted  according  to  the  empirical  Mp  expression  (Fig.  5).  Similar 
behavior  is  observed  for  other  temperatures. 

Interfacial  tensions  for  the  PS/PMMA  system  (Fig.  6)  cover  the  broadest 
range  of  molecular  weight.  The  precision  in  these  measurements  (standard 
deviation  is  ca.  6%)  is  the  lowest  of  the  three  systems  due  to  the  high 
viscosities.  The  data  again  show  apparent  linearity  when  plotted  as  a 
function  of  PS  M 

n 

The  interfacial  tension  data  for  all  three  systems  appear  to  corres- 

-2/3 

pond  well  to  the  empirical  Mn  '  expression  proposed  by  Gaines  and  co- 

Q  “IQ 

workers  at  least  based  upon  the  criteron  of  apparent  linearity  in 
plots  of  y  vs.  Mn”2/^.  A  more  rigorous  estimation  of  the  power  law  of  the 
molecular  weight  dependence  was  obtained  by  performing  non-linear  least 
squares  regression  of  the  data  upon  an  expression  of  form  y  *  ^  -  C£/M*. 
This  analysis  yielded  the  following  values  for  the  exponent  z:  0.54  for 
POMS/PBD  at  25°C;  0.68  for  PS/PBDH  at  1 71 °C;  and  0.90  for  PS/PMMA  at  199°C. 


A  smaller  exponent  for  PDMS/PBD  could  be  explained  by  the  occurence  of 
surface  fractionation  of  the  polydisperse  PDMS  according  to  molecular 
weight.  The  other  materials  are  all  relatively  monodisperse.  Alter¬ 
natively,  the  PDMS/PBD  system  is  closest  to  its  critical  point  and  in 

-0  5 

proximity  to  the  critical  point  an  dependence  of  interfacial  tension 

19  22 

has  been  predicted.  ’  The  PS/PBDH  system  shows  good  correspondence  with 
-2/3 

the  Mn  '  dependence.  A  similar  dependence  for  the  surface  tension  has 

35 

been  explained  using  a  simple  lattice  analysis  that  incorporates  the 

contribution  of  the  end-groups  at  the  interface.  For  these  moderate 

-2/3 

molecular  weights,  the  end-group  effects  are  important  and  an  Mn  depen- 

36 

dence  might  be  expected.  The  PS/PMMA  blends,  on  the  other  hand,  contain 

the  highest  molecular  weight  constituents  and  should  thus  conform  best  to 

the  limit  of  infinite  molecular  weight.  In  this  limit,  the  exponent  z  is 

37 

predicted  to  be  unity. 

The  nonlinear  regression  results,  therefore,  suggest  that  the  exponent 
z  of  the  molecular  weight  dependence  of  polymer-polymer  Interfacial  tension 
increases  as  the  molecular  weight  of  the  constituents  increases.  There  Is 
theoretical  support  for  this  trend,  however  we  should  emphasize  that  the 
precisions  of  the  present  data  are  not  sufficient  to  quantitatively  sub¬ 
stantiate  the  trend.  In  fact  the  quality  of  the  fits  for  the  nonlinear 


regressions  is  not  much  improved  over  the  correspondence  with  the  empirical 
-2/3 

Mn  dependence. 

Comparison  with  Theory 

1 3 

Helfand  and  Tagami  derived  an  expression  for  polymer-polymer  Inter¬ 
facial  tensions  by  performing  a  self-consistent  mean  field  solution  of  a 
segmental  diffusion  equation  across  the  interface.  In  the  limit  of  infinite 
molecular  weight  they  obtained  the  following  closed  form  expression  valid 
for  symmetric  polymers 


1/2 
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where  x  is  the  Flory-Huggins  interaction  parameter 
pQ  the  average  number  monomer  density 
b  the  Kuhn  statistical  segment  length 
k  Boltzmann's  constant;  and 
T  the  temperature  in  K. 

Helfand  and  Sapse14  removed  the  restriction  of  property  symmetry,  and 
obtained 


Y 


kT(V)'/2 

o  2 


Lll i)2 

A  +  6 


(4) 


with  =  pibi  for  species  i  equal  to  A  or  B.  Helfand15  and  Roe16  have 

separately  developed  interfacial  tension  theories  that  are  based  on  lattice 

38 

models  in  the  spirit  of  the  Flory-Huggins  approach.  Experimental  verifi¬ 
cation  of  the  lattice  theories  are  difficult  however,  since  the  parameters 
of  the  lattice  are  unknown  a  priori. 

Comparison  of  our  experimental  data  with  the  predictions  of  the 

Helfand  theories  requires  knowledge  of  pQ,  b  and  x •  pQ  is  obtained  from 

the  experimental  density  data,  while  Kuhn  statistical  segment  lengths  are 

39 

calculated  from  tabulated  characteristic  ratios  . 

There  are  no  literature  data  on  the  Flory-Huggins  interaction  para¬ 
meter  for  the  PBD/PDMS  blend,  and  our  previous  attempts  to  model  this 
system,  employing  the  regular  solution  theory  estimation  of  x  using 
literature  values  of  the  solubility  parameter,  were  not  successful23.  The 
Flory-Huggins  interaction  parameter  and  its  temperature  dependence  for  the 
PS/PBDH  system  was  previously  calculated40  by  small  angle  x-ray  scattering 
analysis  of  homogeneous  poly(styrene-b-hydrogenated  1,2  butadiene)  diblock 


copolymers,  using  the  theory  of  Leibler  .  The  interaction  parameter 
obtained  from  the  analysis  is  given  by 

x  =  xs  +  VT  (5) 

where  xs  =  0.0057,  xH  =  21.0,  and  T  is  the  absolute  temperature. 

This  functional  form  of  the  interaction  parameter  is  consistent  with  the 
notion  that  x  may  be  written  as  the  sum  of  an  enthalpic  and  an  entropic 

30 

contribution.  Refinements  of  the  original  Flory-Huggins  equation  suggest 

that  x  embraces  all  the  changes  in  the  enthalpy  and  the  non-combinational 
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entropy  upon  mixing  .  Non-combinatorial  contributions  to  the  entropy 
of  mixing  arise,  for  example,  from  positive  or  negative  volumes  of  mixing 
and  changes  in  coil  expansion/chain  rigidity  in  the  mixture.  Also  the 
changes  of  order  which  accompany  the  contact  energy  effect  lead  to  a 
non-combinatorial  entropy  of  mixing  contribution.  Expressions  of  the  above 
form  have  already  been  reported  for  polystyrene-polydiene  systems45-48. 

A  comparison  of  the  predictions  of  the  theories  by  Helfand  et  al.  with 
our  experimental  data  appears  in  Fig.  7.  Since  the  theory  explicity 
assumes  infinite  molecular  weight,  we  show  only  the  comparison  for  the 
blend  with  highest  constituent  molecular  weights  (PS  10200/PBDH  4080).  The 
experimental  interaction  parameters  represented  by  (5)  were  used  in  this 
calculation.  The  magnitude  of  the  theoretical  predictions  corresponds  well 
with  the  experimental  data.  The  theoretical  interfacial  tensions  are  found 
to  increase  with  increase  in  temperature,  however,  opposite  to  the  behavior 
of  the  experimental  data. 

Two  of  the  possible  causes  of  the  observed  discrepancy  between  the 
experimental  and  the  predicted  interfacial  tensions  involve  assumptions 
inherent  to  the  calculation.  First  of  all,  we  have  previously  stated  that 


the  theory  explicitly  assumes  infinite  molecular  weight.  Therefore, 
application  of  the  theory  to  our  data  of  relatively  low  molecular  weight 
materials  may  not  be  appropriate.  Secondly,  we  have  assumed  that  the 
interaction  parameter  determined  for  diblock  copolymers  properly  describes 
the  interactions  between  the  corresponding  two  homopolymers  in  an  immiscible 
blend. 

Examination  of  literature  expressions  for  experimental  interaction 
parameters  suggests  that  the  latter  assumption  may  not  be  valid.  Although 
experimental  interaction  parameters  for  binary  homopolymer  blends  and 
diblock  copolymers  both  seem  to  correspond  to  the  relationship  described  in 
Equation  (5),  there  appears  to  be  a  discernible  difference  in  the  behavior 
of  the  temperature  independent  term  for  the  two  types  of  systems.  The 
origin  of  this  term  is  generally  ascribed  to  non-comblnatorlal  entropy 
effects  which  are  not  accounted  for  in  Flory-Huggins  theory,  as  discussed 
earlier  in  this  communication.  There  are  not  many  data  currently  available 
in  the  literature  that  can  permit  detailed  comparison  between  Interaction 
parameters  calculated  from  diblock  copolymers  and  from  the  respective 
homopolymer  blends.  To  our  knowledge  all  xs  values  that  have  been  reported 
for  block  copolymers  have  been  found  to  be  much  smaller  in  magnitude  than 
the  corresponding  values  for  binary  homopolymer  blends. 

A  compilation  of  some  of  the  available  data  is  shown  in  Table  4.  For 

example,  for  the  polystyrene-polyisoprene  system  the  term  varies  from 

-0.9. 10" 3  for  the  binary  blend47  to  +0.1. 10" 3  for  a  similar  diblock  co- 
40  47 

polymer  .  As  pointed  out  by  Rounds  ,  the  parameter  a$  represents  an 
excess  entropy  of  mixing  over  the  combinatorial  entropy  alone.  A  negative 
value  specifies  that  the  mixture  has  greater  excess  entropy  than  the  pure 
components.  From  Table  4,  it  appears  that  this  excess  entropy  is  lower  for 
diblock  copolymers  than  for  binary  polymer  blends.  This  result  may  be 


attributed  to  the  additional  entropy  constraints  intrinsic  to  block  co¬ 
polymers  that  are  associated  with  the  covalent  bond  between  unlike  se¬ 
quences  . 

As  an  alternative  to  using  our  experimental  expression  for  x  from  the 

P(S-b-BDH)  diblocks  in  the  Helfand-Tagami  theory, we  have  applied  the  theory 

directly  to  our  data  to  extract  apparent  x  values  that  would  make  these 

data  meet  the  theoretical  prediction. 

The  apparent  x  values  (Fig.  8)  do  correspond  qualitatively  to  the  form 

of  expression  (5),  but  differ  substantially  from  the  diblock  copolymer  x 

data.  Regression  of  the  apparent  x  parameter  data,  however,  gives  an 

expression  which  is  in  good  agreement  with  data  for  the  other  binary  blends 

(Table  4).  For  the  range  of  temperatures  accessible  to  most  experiments, 

these  apparent  x  values  are  larger  than  those  found  for  PS/PBD,  consistent 

with  the  observation  that  in  diblock  copolymers,  a  hydrogenated  diene 

sequence  (PBDH)  is  less  miscible  with  polystyrene  than  the  original  PBD 
40 

sequence  . 

As  we  have  previously  discussed,  a  thorough  comparison  of  our  results 

with  the  Helfand  theories  is  not  possible,  due  to  the  assumption  of  infinite 

molecular  weight  inherent  to  the  calculation.  The  effect  of  molecular 

weight  on  interfacial  tension  can  be  calculated  however  by  applying  an 
23 

approach  based  upon  an  extension  of  the  square  gradient  approach  of  van 

49  24 

der  Waals  as  modified  by  Cahn  and  Hilliard  . 

Such  an  approach  has  been  used  to  model  the  surface  tensions  of 
50  51-53 

liquids  and  polymer  melts  ,  polymers  at  the  liquid-liquid  interfaces 

54 

in  biphasic  regular  binary  solutions  ,  interfacial  tensions  of  low  and 

high  molecular  weight  liquid  mixtures  and  demixed  polymer  solutions  , 

18  56  57 

and  spinodal  decomposition  in  polymer  blends  *  .  Sanchez  has  shown  that 


the  gradient  theory  is  "in  harmony  with  the  microscopic  theory  of  He! f and 
and  co-workers  although  the  latter  treats  polymer  interfaces  from  a  com¬ 
pletely  different  point  of  view". 

The  basis  of  the  theory  is  the  assumption  that  the  free  energy  per 
unit  volume  in  a  region  of  non-uniform  composition  is  a  function  of  both 
the  local  composition  $  and  the  composition  of  the  immediate  environment. 
The  total  free  energy  for  the  mixture  of  volume  V  is 

G  =  ^  [G0($)  +  < (V4> +....]dV  (g) 

where  Go(<f> )  is  the  free  energy  density  of  a  uniform  system  of  composition 
2 

4> ,  and,  k(v#)  is  the  additional  positive  contribution  to  the  free  energy 
arising  from  the  local  concentration  gradient. 

24 

Following  the  development  by  Cahn-Hilliard,  the  interfacial  tension 
for  the  case  of  an  one-dimensional  composition  change  across  a  flat  Inter¬ 
face  separating  two  coexisting  phases  o  and  B  is  given  as 

*6 

y  =  2  j  (KAg($))^  d<p  (7) 

<Pa,  <|>g  are  the  volume  fractions  of  the  two  coexisting  phases  and  AgU)  is 
the  excess  free  energy  density  of  the  uniform  system  with  respect  to  a 
standard  state  of  an  equilibrium  mixture  of  a  and  b,  i.e. 

ag  ($)  »  aGo($)  -  [nA  Aua  (*e)  +  ng  avb  Ue)]  (8) 
with:  n A  and  ng  the  number  density  of  molecules  of  type  A  and  B;  $  ,  the 

equilibrium  composition  of  either  of  the  coexisting  phases;  and  Aua  and  Aug 

the  changes  In  chemical  potentials  of  A  and  B.  From  the  Flory-Hugglns 
38 

lattice  model  for  the  free  energy,  it  follows  that 


AGo/kT  *  nA  In*  +  ng  ln(l-$)  ♦  x  nA  NA  (1-*) 


(9) 
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Ay./kT  =  In*  +  (1-*)  (1  -  Va  )  +  x  N.  (I-*)2 

Ve 


(10) 


Avig/  kT 
with  * 


=  ln(l-*)  +  *  (1-  Vb)  +  x  Nr  %  * 


NAuA 


nANAuA 


nANAUA  +  nBNBuB 


01) 

.(12) 


Na  and  Ng  are  the  respective  degrees  of  polymerization;  uA  and  ug  the 
specific  monomer  volumes.  The  composition  *  and  *Q  of  the  coexisting 
phase  a  and  e  at  equilibrium  are  found  by  equating  the  chemical  potentials, 
such  that 


(AviA)a  -  Upa)b  and  UwB)a  =  (Aug)6  (13) 

The  coefficient  of  the  square  gradient  form  can  be  derived  using 
linear  response  theory  within  the  framework  of  the  random  phase  approxi¬ 
mation.'*®”^  de  Gennes66  suggested  that  the  coil  remains  Ideal  on  the  scale 
of  one  coil,  even  in  the  case  of  a  dense  mixture  of  interacting  chains. 
Therefore,  an  ideal  single  chain  approximation  can  be  employed  In  the 
calculation  of  the  scattering  function,  S(q),  where  q  Is  the  scattering 
vector.  The  scattering  function  is  related  to  the  volume  fractions  and  the 
chain  length  by60’63 


1  .  1  .  1  & 
S(q)  *uA  N,  fD  (NA,q)  (!-♦)  Ug  Ng  (Ng,q)  uft 

where  fD(N,q)  is  the  Debye  function64,  defined  as 

f0(N,q)  =  2  (U  +  exp(-U)-l) 

U2"” 


with 


U  *  N2q2b2/6  =  q2<r  2>/6  =  q2Rr2 


(14) 


(15) 


(16) 


where  b  is  the  Kuhn  statistical  segment  length,  <r  >  is  the  mean-squared 
end-to-end  distance,  and  RQ  is  the  radius  of  gyration. 


Two  limiting  expressions  for  S(q)  can  be  calculated  for  qRQ  »l  and 
qRQ«l.  The  first  corresponds  to  a  sharp  Interface,  while  the  second 
corresponds  to  a  relatively  diffuse  one.  For  qRr»l,  fn  £  12/Nq2b2  * 


12/q<rQ  >,  and 


1  =  <r02>A  „ 

TqT  VT  tu,  n. 


*UA  NA 


<ro2>B  -  2y 

~7T  r  TT  — 


.2_  2 


whereas,  for  qRg<<l,  f^p-NO-g  <rQ  >  /18),  and  using  the  equation  for  the 


spinodal  curve 


2xj  U) 


iA  uA$  NB  uB  (l-$) 


the  scattering  function  is  given  by 


2  (xs(q)-x)  ,  q2  r<r02»A  +  <r02.B 


18  $uA  Na  <^Ug  Ng 


Therefore  the  square  gradient  coefficient,  <u',  for  the  narrow  Interphase 
is 


<(,)  -„2> 


2  2 
ro  >A  <ro  *8 

^♦u.n.  wm 


and  the  coefficient,  < 


applicable  to  the  broad  interphase  case  is 


2  2 
<  r  > «  <r>„ 

0  A  0  B 

36«TuaNa  3^0-tJUgNg 


Equation  (21)  and  its  equivalent  for  a  symmetric  system  has  been  widely 
used^®’5®’5^’61 ,63  to  model  the  dynamics  of  concentration  fluctuations  in 
binary  polymer  blends  near  the  critical  point,  whereas  equation  (20)  has 
been  used  to  study  micelle  formation  in  homopolymer  copolymer  mixtures. 


A  comparison  between  the  predictions  of  the  square  gradient  theory 
approach  and  our  experimental  interfacial  tensions  for  the  PS  10200/PB0H 
4080  system  is  shown  in  Fig.  9.  The  experimental  values  of  x  given  by  Eq. 

5  were  used  in  the  calculations.  The  predicted  interfacial  tensions 
compare  remarkably  well  with  our  data  for  the  coefficient  of  the  square 
gradient  term,  k,  given  by  Eq.  20,  and  are  within  20%  of  the  experimental 
data  for  the  coefficient  *  given  by  Eq.  21.  Note  that,  no  adjustable 
parameters  were  allowed  in  this  comparison. 

The  theory  was  also  examined  with  respect  to  its  ability  to  predict 
the  molecular  weight  dependence  of  interfacial  tension.  In  Fig.  10  the 
predictions  of  the  theory  are  compared  to  our  PS/PBDH  4080  systems  at  a 
temperature  of  171°C.  At  high  molecular  weights,  using  Eq.  21  for  the 
coefficient  < ,  the  theoretical  curve  corresponds  well  with  the  extrapolated 
empirical  relationship  for  the  experimental  data,  while  use  of  Eq.  20  leads 
to  an  overestimation  of  interfacial  tension  by  ca.  20%.  Adjustable  para¬ 
meters  were  not  allowed  in  this  comparison.  The  theory  does  predict  an 
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apparent  dependence  of  Interfacial  tension  on  Mn  '  *  where  Is  the 
number  average  molecular  weight,  however,  it  deviates  considerably  from  the 
experimental  data  for  low  molecular  weights.  The  predictions  of  the  theory 
are  physically  unrealistic  for  low  molecular  weights,  since  they  erroneously 
indicate  complete  miscibility  (i.e.,  y  =  0)  for  a  PS  molecular  weight  of 
ca.  2400.  During  our  experiments,  however,  two  phases  were  always  present 
under  these  conditions,  and  appreciable  mixing  was  not  observed.  Similarly, 
when  the  theory  was  applied  to  a  PS  2200/polybutadiene  (Mn  =  7800,  M iyM  = 
1.02,  95%  1,2)  blend,  with  interaction  parameter  values  measured  by  Owens40 
from  the  respective  copolymer,  misibility  was  erroneously  predicted  at 
145°C;  two  phases  were  present  however,  during  our  Interfacial  tension 
measurements66. 


The  inability  of  the  theory  to  predict  the  proper  critical  molecular 
weight  does  not  necessarily  indicate  inadequacy  of  the  square  gradient 
approach  but  implies  failure  of  the  Flory-Huggins  free  energy  expression. 
Therefore,  within  the  framework  of  Flory-Huggins  theory,  it  again  appears 
that  the  use  of  x  values  determined  for  diblock  copolymers  is  inappropriate 
in  describing  the  interactions  in  polymer  blends. 

There  may  also  be  deficiencies  associated  with  the  use  of  the  square 
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gradient  theory  itself.  Halperin  and  Pincus  pointed  out  that  because  the 
Cahn-Hilliard  theory  is  a  mean-field  theory,  its  validity  near  the  critical 
point  is  only  qualitative.  On  the  other  hand,  the  theory  assumes  weak 
composition  gradients  that  may  be  realized  only  close  to  the  critical 
region.  Binder*^  suggested  that  for  the  qRg>>]  case,  an  additional 
correction  term  should  be  included  in  the  gradient  term,  which  arises  from 
the  finite  range  of  interactions,  and  is  proportional  to  (a*)  where  i|» 
is  the  range  of  interactions.  ForX<l,  however,  this  correction  term  is 
negligible  as  suggested  by  de  Gennes  .  He  also  argues  that  Eq.  (21)  de¬ 
scribes  well  the  additional  positive  contribution  to  the  free  energy  from 
the  local  concentration  gradients,  even  in  the  case  x>>xc  (xc  is  the  value 
of  the  interaction  parameter  at  the  critical  point),  i.e.,  In  the  strong 
segregation  limit.  Additionally,  polymer  molecules  in  the  interfacial 
region  are  limited  in  their  allowable  configurations,  and  subsequently 
suffer  a  loss  of  configurational  entropy®^.  If  the  interface  is  broad 
relatively  to  the  molecular  size,  these  effects  are  minimal,  and  the  square 
gradient  coefficient  represented  by  Eq.  21,  or  a  slightly  modified  form  of 
Eq.  21  is  appropriate.  In  the  narrow  interface  limit,  however,  both 
entropic  and  energetic  contributions  should  be  included  in  the  gradient 
term.  For  low  molecular  weights,  both  these  contributions  are  important, 
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since  non-local  entropic  contributions  are  comparable  in  magnitude  to  the 
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non-local  enthalpic  contributions  Meier  suggested  that  more 

explicit  Interfacial  constraints  should  be  Included  to  describe  the  loss  of 

configurations  at  the  narrow  interface  limit. 

Even  with  these  uncerta Ini  ties,  however,  we  are  encouraged  by  our 

present  results.  To  our  knowledge,  the  square  gradient  theory  Is  the  first 

-2/3 

thermodynamic  theory  to  predict  an  apparent  Mn  '  dependence  of  inter¬ 
facial  tension,  and  the  first  to  successfully  reproduce  the  temperature 
dependence  of  interfacial  tension. 

Concluding  Remarks 

The  preceding  discussion  points  out  several  problems  associated  with 

the  establishment  of  a  fundamental  understanding  of  Interfacial  tensions  in 

Immiscible  binary  homopolymer  blends.  A  number  of  further  developments, 

both  theoretical  and  experimental,  are  required  to  overcome  these  problems. 

The  most  essential  experimental  requirement  Is  the  knowledge  of  the 

interaction  parameter,  x»  for  the  actual  system  being  studied.  Procedures 

for  the  determination  of  x  parameters  In  binary  homopolymer  blends  have 
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been  discussed  by  Riedl  and  Prud'homme  .  Most  methods  developed  to  date 
require  either  knowledge  of  the  blnodal  curve  of  the  phase  diagram,  or 
accessibility  to  the  miscible  state.  Attainment  of  these  conditions  for 
most  polymer  pairs  restricts  the  materials  to  be  of  rather  low  molecular 
weight.  As  mentioned  earlier,  the  interfacial  tension  measurements  are 
subject  to  this  same  constraint. 

Statistical  thermodynamic  theories,  on  the  other  hand,  have  been 
solved  primarily  for  the  case  of  infinite  molecular  weight.  The  theories 
based  upon  the  square  gradient  framework  can  be  applied  to  low  molecular 
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weight  mixtures,  but  do  not  appropriately  consider  the  loss  of  con¬ 
figurations  available  to  the  chains  in  the  interphase  region.  Thus  in 
order  to  compare  theoretical  predictions  to  accessible  experimental  data, 
one  must  either  solve  the  statistical  thermodynamics  theories  for  finite 
molecular  weight,  or  develop  a  more  rigorous  square  gradient  theory  which 
considers  all  of  the  appropriate  constraints  of  the  chains  in  the  inter¬ 
phase.  These  refined  theories  could  be  compared  directly  to  accessible 
experimental  data  on  polymer-polymer  interfacial  tension. 
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Table  2 

Temperature  dependence  of  Interfacial  Tension  for  PBD  1000/PDMS 

y  =  a  -  b  t 


Polymer  Pair 

a 

b 

(dyn/cm) 

( dyn/cm/ °C) 

PBD  1000/PDMS  770 

1.3409 

1.065.10"2 

PBD  1000/PDMS  1250 

1.5393 

9.381 .10-3 

PBD  1000/PDMS  2000 

2.0647 

7.543.10-3 

PBD  1000/PDMS  3780 

2.8193 

1.190.10*2 
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Table  3 

Temperature  dependence  of  Interfacial  Tension  for  PS/PBDH  4080 

y  s  a  -  b  t 


Polymer  Pair 

a 

b 

(dyn/cm) 

( dyn/cm/ *C) 

PS  2200/PBDH  4080 

3.7139 

1.477.10"2 

PS  4000/PBDH  4080 

3.7468 

1.118.10"2 

PS  10200/PBPH  4080 

3.8329 

8.82. 10'3 
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Literature  Interaction  Parameter  Data 
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Microstructure 

(a  X  1C3)9 
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1,4 

3,4 

(mol /cm3) 

PS/PI3 

- 

90 

10 

-0.9  +  750/T 

47 

P(S-b-I)b 

- 

93 

7 

0.1055+293/Th 

0.0090+25/T 

4C 

P(S-b-I)b 

38 

3 

59 

-0.0937+66/T 

45 

PS/PBD0 

23 

77 

- 

-0.9  +  750/T 

47 

PS/PBDC 

6 

94 

- 

-1 .05  +  537/T 

46 

P(S-b-BD)d 

95 

5 

- 

-0.021  +  25/T 

40 

P(S-b-BD)d 

-0.027  +  28/T 

48 

PS/PBDC 

14 

81 

-0.2024+1 16/T 

18 

P ( S-b-DDH )e 

95 

5 

0.0057  +27/T 

47 

PS/PBDHf 

95 

5 

-1.61+894/T* 

-0.1296+73/TJ 

"Notes 

a.  PS  *  polystyrene;  PI  =  polyisoprene 

b.  P(S-b-I)  =  poly(styrene-b-lsoprene)  diblock  copolymer 

c.  PBD  =  polybutadiene 

d.  P(S-b-BD)  =  poly(styrene-b-butadiene)  diblock  copolymer 

e.  P(S-b-BDH)  =  poly(styrene-b-l ,2  hydrogenated  butadiene)  diblock  copolymer 

f.  PBDH  *  hydrogenated  1,2  polybutadiene  3 

g.  a  =  x  Pm  with  Pu  is  the  average  molar  density  in  rrol/cm3 

h.  calculated  froir.xhe  measured  x  with  pm  =  0.01172  mol/cnr 

i.  calculated  using  the  experimental  densities  and  x 

j.  calculated  using  Helfand-Tagaml  theory  with  experimental  interfacial  tension 
data 
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Fig.  1. 


Fig.  2. 


Fig.  3. 


Fig.  4. 


Fig.  7. 


Fig.  8. 


Fig.  9. 


Quality  of  the  fit  obtained  after  the  application  of  the  algorithm  to  the 
experimental  profile  for  a  PS  drop  in  PBDH.  Solid  line  is  the  theore¬ 
tical  profile  and  the  data  points  denote  the  original  segmented 
experimental  drop  profile. 


Experimental  interfacial  tension  as  a  function  of  temperature  for 
PBD/PDMS  pairs,  open  diamonds:  PBD  1000/PDMS  3780;  filled  squares:  PBO 
1000/PDMS  2000;  filled  diamonds:  PBD  1000/PDMS  1250;  open  squares:  PBD 
1000/PDMS  770 


Experimental  interfacial  tension  between  PDMS  and  PBD  1000  as  a  function 
of  PDMS  M  .  Open  squares,  these  measurements;  filled  squares:  data  taken 
from  Anastasiadis  et  al.  . 


Experimental  interfacial  tension  as  a  function  of  temperature  for  PS/PBDH 
pairs.  Filled  squares:  PS  10200/PBDH  4080;  filled  diamonds:  PS  4000/PBDH 
4080;  open  squares  PS  2200/PBDH  4080 


Fig.  5.  Experimental  Interfacial  tension  between  PS  and  PBDH  4080  as  a 


function  of  PS  Mn  at  1 71 °C- 


Fig.  6.  Experimental  interfacial  tension  between  PS  and  PMMA  10000  as  a  function 


of  PS  M  at  199°C. 
n 


Comparison  of  experimental  interfacial  tension  (points)  for  PS  10200/PBDH 
4080  with  Helfand  theories.  Solid  line:  Helfand-  Tagaml  »  Eq.  4;  Dashed 
line:  Helfand-Sapse  ,  Eq.  5. 


Comparison  of  interaction  parameter  data  of-Owens*0  (diamonds)  with 
values  obtained  from  Helfand-Tagaml  theory  J  (squares)  using  experimental 
interfaclal  tension  data. 


Comparison  of  experimental  Interfaclal  tension  (data  points)  for  PS 
10200/PBDH  4080  with  the  square  gradient  theory.  The  square  gradient 
coefficient  is  given  by  Eq.  20  (solid  line)  and  Eq.  21  (dotted  line). 


Fig.  10. 


Comparison  of  experimental  Interfaclal  tension  (data  points)  for  PS/PBDH 
4080  at  171°C  with  the  square  gradient  theory.  The  square  gradient 
coefficient  Is  given  by  Eq.  20  (solid  line)  and  Eq.  21  (dotted  line). 
The  dashed  line  is  the  linear  fit  of  the  data. 
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